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Recent  experiments  suggest  that  the  isomerization  reaction  CH30  =*  CH2OH  may  compete  with  the  hydrogen  dissocia¬ 
tion  reaction  CH3O  -*  CH^O  +  II.  We  report  correlated  many-body  calculations  of  barrier  hqjghts  for  these  reactions  and 
other  decompositions  of  CH3O  and  CH2OH.  The  barriers  for  isomerization  and  dissociation  are  similar,  while  the  reaction 
CH3O  —  HCO  +  II 2  has  a  higher  energy  barrier. 


I .  Introduction 

The  methoxy  radical  CH3O  is  thought  to  play  an 
important  role  in  hydrocarbon  oxidation  processes 
relevant  to  combustion  and  photochemical  air  pollu¬ 
tion  [1] .  Recently,  Radford  [2]  has  suggested  that  the 
unimolecular  isomerization  of  methoxy  to  form  the 
hydroxymethylene  radical, 

CH30-CH20H,  (1) 

followed  by  the  fast  reaction 

ch2oh  +  o2  -ch2o  +  ho2,  (2) 

may  provide  an  alternate  route  for  the  oxidation  of 
methyl  radicals.  Subsequently,  Batt  et  al.  [3]  pointed 
out  that  although  the  isomerization  reaction  should  be 
competitive  with  the  hydrogen-dissociation  reaction 
of  the  methoxy  radical, 

CH30~>>CH20  +  H,  (3) 

neither  reaction  should  occur  quickly  enough  to  com- 

*  This  work  was  supported,  in  part,  by  the  Army  Research  Of¬ 
fice  under  contract  No.  DAAG29-80-C-0105  to  the  Universi¬ 
ty  of  Florida. 


pete  with  other  faster  reactions  of  the  radical  that  are 
important  in  atmospheric  chemistry.  However,  at  the 
higher  temperatures  encountered  in  combustion  pro¬ 
cesses,  the  unimolecular  reactions  may  be  more  signifi¬ 
cant. 

Here  we  report  computed  classical  barriers  for  sever¬ 
al  unimolecular  reactions  of  CH30  and  CH2OH.  In  ad¬ 
dition  to  (1)  and  (3),  we  also  discuss  the  hydrogen-dis¬ 
sociation  reactions  of  hydroxymethylene, 

CH20H^CH20  +  H,  (4) 

and  the  hydrogen  molecule  elimination  reaction  of 
methoxy, 

CH30-*CH0  +  H2.  (5) 

2.  Outline  of  the  calculations 

For  each  of  the  open-shell  species  in  this  study  an 
unrestricted  Hartree-Fock  (UHF)  wavefunction  is 
chosen  for  the  ground  electronic  reference  state.  The 
structural  parameters  for  the  reactants  and  the  various 
transition  states  are  determined  using  the  analytic  gra¬ 
dient  based  GRADSCF  codes  [4]  and  a  6-3 1 G**  basis 
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(5) .  The  many-body  perturbation  theory  (MBPT)  cal¬ 
culations  use  double  zeta-quality  basis  sets  [6]  aug¬ 
mented  with  polarization  functions  on  all  atoms  [6] . 
Since  self-consistent-field  calculations  provide  poor  es¬ 
timates  of  energy  differences,  the  electronic  energy  is 
determined  using  MBPT  [7| ,  In  particular,  fourth-order 
MBPT  calculations  that  include  all  single-,  double-,  and 
quadruple  excitation  diagrams  [SDQ  MBPT(4)|  are 
reported. 


3.  Equilibrium  configurations  for  methoxy  and 
hydroxymethylene 

The  structural  parameters  determined  for  the  meth¬ 
oxy  and  hydroxymethylene  radicals  are  shown  in  fig. 

1.  Our  SCF  calculations  predict  structural  parameters 
for  CH30  that  agree  closely  with  those  of  Yarkony  et 
al.  [8] .  Each  of  the  SCF  calculations,  however,  dis¬ 
agrees  somewhat  with  the  results  of  our  previous  cor¬ 
related  D  MBPT(4)  calculation  for  this  radical  (9] .  In 
that  study.  Bent  et  al.  showed  that  the  C3v  symmetry 
axis  is  broken  by  a  reduction  in  the  in-plane  OCH  bond 
angle,  while,  to  the  contrary,  the  SCF  calculations  pre¬ 
dict  that  the  symmetry  changes  via  an  increase  of  the 
in-plane  OCH  bond  angle.  A  correlated  calculation  of 
the  CH-jO,  CH3OH.  CH20.  and  HCO  structures  has 
also  been  published  [10| .  Since  the  structur  „  of  the 
radicals  and  transition  states  described  in  this  work  are 
determined  using  UHF  calculations,  for  consistency, 
all  energy  differences  for  the  methoxy  reactions  are 
computed  at  the  structural  parameters  determined  in 
the  UHF  calculations. 

The  electronic  energies  for  the  two  radicals,  includ¬ 
ing  both  the  UHF  and  SDQ  MBPT(4)  results,  are  re¬ 
ported  in  table  1.  Interestingly,  the  UHF  calculations 
suggest  that  the  methoxy  radical  CH30  is  more  stable 


IIOM* 


CHjOH  (eq) 


Fig.  1.  UHF  predicted  equilibrium  structures  for  CH3O  and 
CHqOH  (bond  lengths  in  bohr). 
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Table  1 

Electronic  energy  predictions  (hartree) 


Structure 

UHF 

SDO  MBPTI4) 

CHjO,  equilibrium 

114.4484 

1 14.7621 

CHjOII.  equilibrium 

114.4419 

114.7684 

1, 2-shift  SPa) 

-114.3668 

1 14.7054 

CH30-CH20  +  H  SP 

114.3845 

114.7060 

CHjOH  -*  CHjO  +  II  SP 

1  14.3679 

114.6923 

a*  SP  =  saddle  point. 

than  the  hydroxymethylene  radical  CH2OH.  The  ener¬ 
gy  difference  equals  4.07  kcal/mole.  On  the  other  hand, 
the  correlated  SDQ  MBPT(4)  calculations  predict  that 
the  hydroxymethylene  radical  is  more  stable,  by  3.90 
kcal/mole.  Experimental  studies  [1 1.12]  indicate  that 
hydroxymethylene  is  5- 10  kcal/mole  more  stable  than 
methoxy.  although  a  more  recent  experimental  [13] 
determination  of  the  heat  of  formation  of  the  methoxy 
radical  would  reduce  this  range  by  3  kcal/mole.  Our 
result  of  3.90  is  in  good  agreement  with  this  lower 
<2—  7  kcal/mole)  estimate. 


4.  Isomerization  and  dissociation  reactions 

The  isomerization  reaction  that  converts  methoxy 
to  hydroxymethylene  belongs  to  the  reaction  class  iden¬ 
tified  as  1 .2 -hydrogen  shifts.  A  recent  review  by  Schaefer 
[14]  discusses  this  class  of  reactions.  The  transition  state 
identified  for  the  isomerization  reaction  is  shown  in 
fig.  2.  Such  isomerizations  are  often  ignored  in  chemical 
kinetic  models,  since  it  is  assumed  that  the  potential 
energy  barrier  for  the  reaction  is  large.  As  the  results 
in  table  1  show,  however,  the  electronic  energy  for  this 
transition  state  does  not  differ  significantly  from  those 
values  calculated  for  the  two  hydrogen-dissociation  re¬ 
actions. 

In  particular,  as  summarized  in  table  2.  the  classical 
barrier  for  the  isomerization,  reaction  (1).  is  35.6  kcal/ 
mole  (39.5  kcal/mole  for  the  reverse  reaction)  with  the 
competing  H-dissociation,  reaction  (3),  having  a  35.2 
kcal/mole  barrier.  The  UHF  transition  state  structures 
arc  shown  in  fig.  2.  Since  the  1 .2  shift  transition  state 
has  a  more  compact  structure,  one  would  expect  that 
the  zero-point  vibrational  energy  differences  would  re- 
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H-elimmotton,  CHjO  H-elimiAolion,  CHjOH 


l-'ig.  2.  UHL'  predictions  of  transition  states  for  C1I30  under¬ 
going  1,2-H  shift.  II  elimination,  Hj  elimination,  and  for  H 
elimination  from  CH2OH  (bond  lengths  in  bohr). 

suit  in  a  lower  0  K  activation  barrier  for  (3)  than  for 
( I ).  Nevertheless,  the  classical  barriers  are  so  nearly 
equal  that  it  is  expected  that  both  reactions  will  occur 
under  most  conditions. 

In  addition  to  the  above,  the  H-dissociation  reaction, 
reaction  (4),  is  also  of  interest.  The  competing  reaction 

CH2OH  -»HCOH  +  H  (6) 

is  ruled  out  since  its  energy  barrier  is  greater  than  85 
kcal/mole  1 15-17] . 

The  saddle-point  structure  for  (4)  is  shown  in  fig.  2. 
SDQ  MBPT(4)  calculations  at  the  UHF  determined 
geometries  suggests  a  barrier  of  47.8  kcal/mole  for  (4). 
Consequently,  it  would  appear  that  the  preferred  uni- 
molecular  reaction  of  CH2OH  is  the  isomerization  to 

ch3o. 


Table  2 

Reaction  barriers  (kcal/mole) 

Reaction  Barrier  height 

CHjO  -  CHjO  +  H  35.2 

CH30~0!2OH  35.6 

CHs0H-CH30  39.5 

CHjOH-CHjO  +  H  47.8 

a*  All  barrier  heights  are  classical.  No  attempt  has  been  made 
to  correct  for  differences  in  the  zero-point  vibrational  ener¬ 
gies. 


5.  The  H2  -elimination  reaction  of  methoxy 

The  reaction  of  methoxy  to  produce  the  formyl 
radical  and  hydrogen  molecule  is  similar  to  the  molec¬ 
ular  products  reaction  channel  of  formaldehyde  (16 
22] .  Like  formaldehyde,  the  methoxy  radical  has  a 
hydrogen-dissociation  channel,  an  isomerization  channel, 
and  an  H2 -elimination  channel.  Recent  MINDO/3  re¬ 
sults  suggest  that  the  dominant  dissociation  of  the  ra¬ 
dical  is  the  H2 -elimination  reaction  [23] .  Furthermore, 
those  results  suggest  that  the  isomerization  reaction  and 
the  H2-elimination  reaction  have  nearly  the  same  activa¬ 
tion  energy. 

We  have  made  numerous  efforts  to  locate  the  H2- 
elimination  reaction  saddle  point  using  the  GRADSCF 
computer  codes.  The  structure  shown  in  fig.  2  corre¬ 
sponds  to  a  structure  for  which  the  norm  of  the  gradient 
satisfies  a  stringent  convergence  criterion.  Proper  iden¬ 
tification  of  a  saddle  point,  however,  requires  the  calcu¬ 
lation  of  the  vibrational  frequencies  for  the  structure. 

If  one  of  the  vibrational  frequencies  is  imaginary,  then 
the  structure  corresponds  to  a  saddle  point.  The  struc¬ 
ture  shown  in  fig.  2  fails  this  test,  because  the  vibra¬ 
tional  analysis  yields  two  imaginary  frequencies.  Addi¬ 
tional  study  of  the  hypersurface  in  this  region,  how¬ 
ever,  suggests  that  the  second  imaginary  frequency  is 
due  to  an  artifact  in  the  SCF  results  for  this  structure. 

In  particular,  our  UHF  wavefunction  breaks  the  Cs  sym¬ 
metry  assumed  for  the  problem.  The  symmetry  break¬ 
ing  may  be  viewed  as  due  to  the  existence  of  a  large- 
amplitude  wagging  motion  characterized  by  a  double¬ 
minimum  potential  energy  function.  This  vibration  may 
be  viewed  as  a  rocking  of  the  H2  fragment  perpendicu¬ 
lar  to  the  plane  defined  by  the  formyl  radical  fragment. 

To  test  this  view,  we  extended  the  “transition  state"’ 
along  the  apparent  reaction  coordinate,  and  searched 
for  a  minimum  on  the  hypersurface.  This  calculation 
leads  to  a  minimum  that  corresponds  to  the  products 
H2  and  HCO.  Distortion  of  the  structure  along  the 
normal  mode  of  the  second  imaginary  frequency  fol¬ 
lowed  by  steepest-descent  minimization  appears  to  lead 
to  no  products.  Instead,  the  search  routine  wanders  in 
hyperspace  near  the  transition  state  structure.  After 
twelve  evaluations,  the  structure  remains  almost  equiva¬ 
lent  to  that  located  by  the  transition  state  search  rou¬ 
tine. 

Finally,  to  test  whether  the  symmetry  breaking  was 
an  artifact  of  the  UHF  calculation,  a  transition-state 
search  was  carried  out  using  an  analytic  gradient  based 
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restricted  Hartree-Fock  calculation.  This  calculation 
also  converged,  in  forty-eight  function  evaluations,  to 
a  structure  with  no  plane  of  symmetry.  The  bond 
lengths  and  bond  angles  predicted  in  that  calculation 
are  close  to  those  predicted  by  the  UHF  calculations. 
Furthermore,  when  force  constants  for  the  structure 
are  calculated  at  the  RHF  level,  we  still  obtain  two 
imaginary  frequencies.  Thus  there  is  no  substantial  dif¬ 
ference  between  the  UHF  and  RHF  predictions. 

In  sum,  our  calculations  to  date  demonstrate  that 
the  transition  state  does  lead  to  the  hydrogen  molecule 
product.  The  calculations  do  not  conclusively  demon¬ 
strate  the  existence  of  a  double-minimum  potential  well 
about  the  transition  state.  To  determine  whether  this 
potential  is,  as  we  believe,  an  artifact  of  the  SCF  level 
calculations  requires  a  re-examination  of  this  transition- 
state  region  using  a  theoretical  method  that  includes 
correlation  effects.  We  anticipate  reporting  such  calcu¬ 
lations  soon. 

The  electronic  energy  including  the  correlation  en¬ 
ergy  for  this  transition  state  is  -1 14.6732  au,  and  the 
reaction  barrier,  55.8  kcal/mole,  can  be  compared  with 
the  other  reaction  barriers  in  table  2.  The  harrier  for  the 
H2 -elimination  reaction  exceeds  the  barriers  determined 
for  the  other  methoxy  radical  reactions  by  more  than 
20  kcal/mole. 


6.  Discussion  and  conclusion 

The  energy  barrier  results  determined  by  these  cal¬ 
culations  and  summarized  in  table  2  support  Radford's 
contention  that  the  hydrogen-dissociation  and  isomeri¬ 
zation  reactions  of  methoxy  radical  could  compete 
favorably  under  appropriate  conditions.  In  addition, 
these  results  suggest  that  for  hydroxymethylene,  the 
isomerization  reaction  should  be  faster  than  the  hydro¬ 
gen-dissociation  reaction.  Finally,  our  preliminary  re¬ 
sults  indicate  that  the  Hj-elimination  reaction  of  the 
methoxy  radical  has  a  much  larger  energy  barrier  than 
either  of  the  other  methoxy  reactions  considered  here. 
In  all  cases,  the  calculations  do  not  consider  the  some¬ 
times  important  effect  of  the  fourth-order  triplet-exci¬ 
tation  diagrams  [22,24,25]  and  are  limited  to  UHF 
determined  structures.  For  example,  the  triple-excita¬ 
tion  contribution  in  H2CO  -*  H2  +  CO  reduced  the  bar¬ 
rier  by  ®=3  kcal/mole  [22]  and  it  is  possible  to  obtain 
similar  changes  by  using  better,  correlated  structures. 
The  spin  contamination  present  in  UHF  calculations 


can  sometimes  cause  additional  problems.  For  the  H2 
elimination  pathway,  the  UHF  multiplicity  is  2.013 
and  SDQ  MBPT(4)  improves  this.  We  expect  that  a 
correlated  determination  of  the  (5)  transition  state 
should  eliminate  any  possible  artifacts  with  the  SCF 
transition  state. 
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